Multiphoton ionization and stabilization of helium in superintense xuv fields 
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Multiphoton ionization of helium is investigated in the superintense field regime, with particular emphasis on 
the role of the electron-electron interaction in the ionization and stabilization dynamics. To accomplish this, we 
solve ab initio the time-dependent Schrodinger equation with the full electron-electron interaction included. By 
comparing the ionization yields obtained from the full calculations with corresponding results of an independent- 
electron model, we come to the somewhat counterintuitive conclusion that the single-particle picture breaks 
down at superstrong field strengths. We explain this finding from the perspective of the so-called Kramers- 
Henneberger frame, the reference frame of a free (classical) electron moving in the field. The breakdown 
is tied to the fact that shake-up and shake-off processes cannot be properly accounted for in commonly used 
independent-electron models. In addition, we see evidence of a change from the multiphoton to the shake-off 
ionization regime in the energy distributions of the electrons. From the angular distribution it is apparent that 
correlation is an important factor even in this regime. 

PACS numbers: 32.80.Fb, 32.80.Rm 
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I. INTRODUCTION 

More than 20 years ago, theoretical studies of atomic hy- 
drogen in ultraintense, high-frequency laser fields produced 
an unexpected result 11H9D: when increasing the intensity of 
the laser pulse to such a degree that the applied forces domi- 
nate over the Coulomb attraction between the nucleus and the 
electron, the ionization probability does not increase accord- 
ingly, but rather stabilize, or start subsiding. This counterin- 
tuitive phenomenon was dubbed atomic stabilization, and was 
subject to much research in the following decade. The discus- 
sions, controversies and conclusions are available in a num- 
ber of review articles, see e.g. [ilOl4l2ll and references therein. 
It has also been pointed out that atomic stabilization has a 
classical counterpart ifTsl [Till . (See also ifioll and references 
therein.) 

At the start of the nineties, the laser technology required 
to experimentally observe the stabilization effect in tightly 
bound systems was not available. For example, in order to 
measure stabilization in atomic hydrogen, photon energies ex- 
ceeding 13.6 eV, the binding energy of the atom, and intensi- 
ties on the order of 1 W/cm^ or more are required [ fisl [l6l . 
Grobe and Eberly iflTll demonstrated that stabilization could 
occur in at moderate intensities (^ lO^-' W/cm^) and pho- 
ton energies (^ 2 eV), and Wei et al. |18] suggested an ex- 
periment, in which a laser, of realistic frequency and inten- 
sity, could possibly stabilize the unstable He^ ion. How- 
ever, at present, the only experimental confirmations of sta- 
bilization are from studies of low-lying Rydberg states (HH- 
l22fl . With recent advances in free-electron laser (PEL) tech- 
nology, extremely high peak intensities have been achieved, 
with wavelengths ranging from vacuum ultraviolet to soft x- 
rays ll23i l24l] . and even higher intensities are expected to be 
delivered in the near future ll25ll . Thus, laser technology is ap- 
proaching the regime needed for observing atomic stabiliza- 
tion in ground state (neutral) atomic systems. 

Although atomic stabilization has been studied extensively 
during the last two decades, studies of stabilization in systems 
containing two electrons are still scarce lUoi |26| and most of- 



ten assessed with simplified physical models of reduced di- 
mensionality lfT7ll27l - l29ll . Including a second electron adds 
a new dimension to the problem, manifested through the 
electron-electron repulsion. The studies mentioned above re- 
vealed that the electron-electron interaction suppresses atomic 
stabilization in two-electron systems. Although ab initio cal- 
culations of helium have previously been performed at fairly 
high intensities in the xuv regime [30, 31], only recently was 
such endeavors extended into the stabilization regime (HH, 
confirming the detrimental effect of the electron-electron in- 
teraction on stabilization. However, it was shown that the ef- 
fect is markedly less than predicted in models of reduced di- 
mensionality. 

In this work, we revisit the problem of multiphoton ion- 
ization of helium in superintense high-frequency fields. In 
continuation of the work of Birkeland et al. Il32ll we look 
more closely into the strong-field ionization dynamics of the 
atom, with paiticular emphasis on atomic stabilization, con- 
sideiing laser pulses of various central frequencies and dura- 
tions. A comparison of the ionization yields obtained from the 
ab initio calculations, including correlation, with correspond- 
ing results obtained from an independent-electron model re- 
veals that the validity of the latter breaks down at strong fields. 
An analysis of the system equations in the so-called Kramers- 
Henneberger frame [33-36] shows that the electron-electron 
interaction plays a decisive role in this limit. We further show 
that this is manifested in the energy and angular distributions 
of the ejected electrons. 

Atomic units, where nie, h, and e are scaled to unity, are 
used throughout unless stated otherwise. 



II. METHODS 
A. Ab initio calculations 

We obtain the ionization probability of ground state helium 
in extreme laser fields from first principles, i.e., by solving 
(numerically) the full time-dependent Schrodinger equation 
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(TDSE). Formulating the problem in the velocity gauge, the 
Hamiltonian assumes the form 

A sine-squared carrier envelope was chosen for the laser in- 
teraction, 

A^(f) =Aosin2(y)cos(a)f). (2) 

Here Aq = Eq/ co, Eq is the amplitude of the electric field, co 
is the laser frequency, and T is the total pulse duration. The 
semiclassical treatment of the field is a valid approach due to 
the enormous photon flux of superintense lasers. The pulse 
fulfills the constraint of a physical pulse 1,37.] . 

rE{t)dt^Q. (3) 
Jo 

Propagation and analysis of the wavefunction is performed 
with the PyProp framework f3§\, a Python/C-H- software 
package for solving the TDSE. 

The wavefunction is expanded in a B-spline basis 
for each of the radial components, and a coupled spherical 
harmonic basis for the angular components, 

•^(r^r^,^)- L^0i^^^f^(ii.,ii2), (4) 

where k ~ {L,M,li,l2} is a combined index for the angular 
indices. The coupled spherical harmonic basis functions, 

^^f;^(ni,i22) = J^{hl2mM-m\LM)Yl^{ai)Yl^-'"{Cl2) 

m 

(5) 

are obtained by linearly combining products of ordinary 
spherical harmonics, weighted by Clebsch-Gordan coeffi- 
cients |41]. 

As the B-spline basis functions are not orthogonal, an over- 
lap matrix Sij ~ f Bi{r)B j{r)dr is introduced for each elec- 
tronic coordinate. From these the total overlap matrix is found 
for every angular momentum component by taking the Kro- 
necker product S = I^; S i (g) S2, where \ denotes the identity 
matrix and k is the angular index. The resulting TDSE may 
then be written as 

/S^c(f ) = H(f )c(f) (6) 

on matrix form. 

We solve the TDSE using a scheme based on the first-order 
approximation to the matrix exponential 

exp(-/AfS"'H) ==I-/AfS"'H + C>(Af2). (7) 

A direct application of this formula is not desirable due to nu- 
merical instabilities. Instead, we combine one half-step for- 
ward in time, 

c(f+Af/2)= ('l-^S-iH')c(f), (8) 



with one half-step backward in time, 

c(f + Af/2)=. (^I + ^S-'H^c(f+Af), (9) 

to obtain the unconditionally stable Cayley-Hamilton form of 
the time propagator, 

(^S + ^H)c(? + Af)= (^S-|^H^c(0. (10) 

This linear system of equations is too large to be solved di- 
rectly, hence we use an iterative method. Since the matrix 
(S + ^H) is not Hermitian, our choice is the generalized min- 
imum residual method (GMRES), a Krylov subspace method 
which combines Arnoldi iterations with a least-squares prob- 
lem in the projected space f42!,'43'|. In the GMRES algorithm 
the error in the least-squares residuals is controlled by the di- 
mension of the Krylov subspace, which can be increased until 
the desired precision is obtained. 

B. Calculating ionization 

In this work we compute the ionization probability resolved 
in direction and energy. We also do a series of smaller sim- 
ulations, calculating only the total ionization probabilities. 
Separating the single and double ionization is achieved by 
projection onto double continuum states. In order to obtain 
these continuum states exactly one needs to solve a scattering 
problem for the full two-particle system. As this is computa- 
tionally cumbersome, an approximation using single-particle 
states is adopted instead. It can be described as follows. In the 
case of double ionization, when both electrons are far from the 
nucleus, a product of continuum He+ (Z = 2) states is used. 
For single ionization, when one electron is close to the nucleus 
and the other far away a product of bound He+ and continuum 
H(Z=l)is used iH. 

The single-electron states are not orthogonal to the bound 
states of the two-electron system, which may become popu- 
lated during the action of the pulse. Therefore, the projection 
of the final wavefunction on the doubly bound states are re- 
moved before further analysis is conducted. Moreover, as the 
electron-electron correlation is neglected in the double contin- 
uum states, the system must be propagated after the pulse for 
all quantities to converge ll45ll . 

On the other hand, when only calculating the total ioniza- 
tion, a small radial box is sufficient. It is no longer necessary 
to propagate the system after the pulse, in order to minimize 
the interaction term, nor to project onto continuum He+ states. 
An absorbing potential is applied at the box boundary in or- 
der to absorb the emitted electrons and to minimize reflection. 
When coupled with an absorbing potential, we find that only 
about a third of the radius needed to resolve the differential 
probabilities is necessary. The total ionization probability is 
simply the complement of the probability of being in one of 
the bound states. 

To find the bound states, we use the implicitly restarted 
Arnoldi method |46]. This is a version of the Arnoldi method 
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for finding eigenpairs that refines the Krylov subspace basis in 
order to find the wanted eigenvectors and eigenvalues. As the 
Amoldi method tends to find the largest eigenvalues, we also 
use shifted inverse iterations, which let us find the eigenvalues 
near a given value. 

Further details on the discretization, the time integration, 
and the analysis were presented in a recent communica- 
tion IH. 

C. Independent-electron model 

In order to gauge the importance of the electron-electron 
interaction, we repeat the calculations using an independent- 
electron (IE) model [48]. The total wavefunction is approxi- 
mated as a product of two single-electron wavefunctions 

1'(ri,r2) = v^sAE(ri)rHe+(r2)- (H) 

The subscript SAE refers to the single-active electron approx- 
imation. This is a common approximation for many-electron 
problems, which focuses on one electron at the time. Any de- 
pendence on the rest of the electrons is included in a common 
potential that is constant with regards to the other electron po- 
sitions. To find the first electron wavefunction (/SAE- we apply 
a pseudo potential, which includes the shielding of the nucleus 
caused by the other electron ll49ll . 

Z + aye-"v + a^re-''*'- + ase-"(''' 

V(r) = • (12) 

r 

For helium the effective charge Z = 1 and the coefficients 
fli ^ 1.231, fl2 = 0.662, fl3 1.325, 04 = 1.236, as = 0.231, 
and = 0.480 were adopted. The other electron moves in a 
He+ potential, and it is therefore an accurate model for the 
singly ionized atom. The IE model reproduces the correct 
ground state energies and single and double ionization thresh- 
olds, and decently represents the excited states. As the name 
of the model suggests, the electrons do not interact with each 
other, beyond what is included in the shielded nuclear poten- 
tial. That makes this a three-dimensional, rather than a six- 
dimensional problem, and it can be calculated with relative 
ease on an ordinary computer As a consequence of working 
with independent particles, the total (single + double) ioniza- 
tion probability becomes 

^total = 1 ^^SAE^He+ (^■^) 

where Pg^g and are the probability of the SAE and the 
He+ electron, respectively, being in a bound state. The prob- 
ability for double ionization is obtained from the product 

pion _ pion pion t'\A\ 
Mouble — ^SAE^He+ ^^^> 

and the single ionization probability is 

pion _ pb pion .pion pb /ic\ 
^single — ^SAE^He+ ^SAE^He+ ' ^^-^> 

where Psae ^He"+ ionization probabilities of the 

SAE and the He+ electrons. 
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FIG. 1 : (Color online) Ionization probabilities plotted as functions of 
the electronic displacement (Eq/co^) for the frequencies o) = 4a. u. 
(left), CO = 5a.u. (middle) and ft) = lOa.u. (right). In each panel the 
pulses are of 3, 6, 12, and 24 cycles duration from bottom to top. 
Dashed lines: IE calculations. 



D. Convergence of the calculations 

When doing the largest calculations, the radial domain typ- 
ically extends to 80 a. u., although the double of this was em- 
ployed to test the convergence. A 5th order B-spline basis 
of 185 splines is used, distributed exponentially near the nu- 
cleus, and linearly further away. Up to 300 splines were used 
for convergence test purposes. Regarding the angular basis of 
coupled spherical harmonics, / < 5 and L<6 was found to be 
sufficient. Note that the system retains cylindrical symmetry 
in the presence of the z polarized laser field. Therefore, the M 
quantum number is set to throughout. Based on the calcula- 
tions with a larger basis, the error is estimated to be less than 
1% in the ionization probabilities. 

For the smaller calculations, intended to provide only the 
total ionization probability, we use a small radial box of 30 a.u. 
and 80-100 B splines of order 7, distributed linearly. Note 
that we have only a third the box size, but half the number 
of B splines. In these calculations the angular basis went up 
to Z = 7 and L = 6. The small box made it possible to go to 
higher intensities and pulse lengths than for the large box. The 
error in the ionization probability is gauged to be less than 5% 
when Eq/ox^ > 1 a.u., and less than 2% for lower intensities. 



III. RESULTS 

A. Total ionization probabilities 

Figure [T] shows the total (single + double) ionization prob- 
ability versus oCq for three different laser frequencies, co — 4- 
(left panel), 5 (intermediate panel) and 10 a.u. (right panel), 
and for four different pulse durations, 3, 6, 12 and 24 cy- 
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FIG. 2: (Color online) Ionization probabilities for a constant pulse 
duration of Ina.u. The lines correspond to laser frequencies of 
ft) = 4,6, 8, 10a. u., from top to bottom, or equivalently pulse lengths 
of 4, 6, 8, 10 cycles. The dashed lines are the IE calculations. 



cles (from bottom to top). Notice that on the abscissae, the 
domains are given in Oq = Eo/(0^, instead of intensity or 
peak electric field strength. Here Oo represents the displace- 
ment amplitude of a free classical electron in the oscillating 
field |4J. This scaling allows us to easily compare the re- 
sults obtained with different laser frequencies. In most of 
the considered cases, the ionization probability increases with 
Oo up to some point, where it attains a maximum before it 
starts to decline, i.e., we are entering the so-called stabiliza- 
tion regime. When stabilization occurs, the ionization peak 
(corresponding to the 'death valley' lHoll ') is typically situated 
between Oo = 0.6 and 0.7 a. u., independent of laser frequency 
and pulse duration. For very short pulse durations, e.g. the 3- 
cycle pulse of © = 10 a. u., we observe a knee in the function, 
rather than a peak at the stabilization point. This is probably 
due to the relatively large bandwidth of these short pulses and 
the averaging this leads to. For long pulse durations, e.g. the 
24-cycle pulse of co — 4a.u., the atom is almost fully ionized, 
and the stabilizing effect turns out to be weak. The dashed 
lines in the figure are the results of the independent-electron 
model. They show good agreement with the full calculations 
for weak fields (Oo < 0.5 a. u.), and for long pulses, but oth- 
erwise tend to overestimate the stabilization. As a matter 
of fact, the results show that the electron-electron interaction 
suppresses stabilization in all cases. We will return to the rea- 
son for this later 

Figure |2] shows the ionization probability as a function of 
OCo for a pulse of constant duration T = 2;ra.u., but for vary- 
ing frequencies, co — 4, 6, 8, and 10 a.u. The corresponding 
results of the IE model are shown in dashed lines. One im- 
mediately perceives that for higher frequencies, the atom sta- 
bilizes at lower ionization probabilities, in accordance with 
the results in Fig. [T] Note that in the limit of weak fields 
single ionization is by far the dominating ionization channel. 



Thus, from first order perturbation theory (XqT. Now, 

since the pulse duration is kept fixed in Fig.|2](as opposed to 
Fig. [T]), this explains why the results of the calculations with 
different frequencies almost coincide at smaller fields. The 
figure also demonstrates the fact that the discrepancy between 
the IE model (dashed lines) and the full calculations (solid 
lines) increases with the intensity. Furthermore, the stabiliz- 
ing effect turns out to be very weak in the fully correlated 
system. Whereas the full ab initio calculations give ioniza- 
tion probabilities that level off (low frequencies) or increase 
(high frequencies) for high intensities, the IE model returns 
probabilities that are noticably lower As the intensity grows, 
so do the discrepancy. As such, the simplified model tends 
to always underestimate the ionization probability, with the 
consequence that the stabilization effect is overestimated. 

B. The role of electronic correlation 

Figures[T]and|2]clearIy demonstrate that the validity of the 
independent-electron model (fTTl i breaks down in the superin- 
tense field regime. This may appear counterintuitive, as one 
might well expect the opposite to happen, i.e., that the im- 
portance of the electron-electron interaction should be negli- 
gible in the presence of a strong external perturbation. The 
reason why the electron-electron interaction in fact becomes 
more important in this limit can be understood by analyz- 
ing the dynamics in the so-called Kramers-Henneberger (KH) 
frame |33|-|36']. the rest frame of a classical free electron in the 
laser field. In this frame the Hamiltonian, Eq. ([T]i, is cast into 
the form, 

HKH^Y^(^+VKH[r. + a{t)])+^^— (16) 
,ti V 2 / |ri-r2| 

where 

is the Kramers-Henneberger potential, and 

a{t) = f A,[t')dt'z (18) 
Jo 

represents the position relative to the laboratory frame of a 
classical free electron in the electric field E^{t) = —dA,/dt. 
One characteristic feature of the KH frame is that the dipole 
interaction terms enter into the electron-nucleus Coulomb po- 
tentials [c.f. Eq. ([TtIiI. which in turn become time-dependent 
and modified by the external field. Note also that the electron- 
electron interaction term is left unaffected by the HK trans- 
formation. Assuming for the moment that the Hamiltonian 
is periodic in time, i.e., neglecting the pulse profile, the KH 
potentials, Eq. ( fTTI l. are expanded in a Fourier series as 

yKH[r, + a(f)] -^y„(ao,r,)e-'"™, (19) 

n 

with 

K(ao,r,) = ^ f e-'"'^'VKH[ri + a{t)]dt. (20) 
1 Jo 
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Inserting the expansion ( fT9] l into the TDSE and apply- 
ing high-frequency Floquet theory, Gavrila et al. LIO^ .50i-i52ll 
showed that the « = component in Eq. (|20] | plays an in- 
creasingly important role in the dynamics at higher value of 
Oq. Furthermore, in the limit of superintense fields (c«() 3> 
1), F0rre et al. [15] showed that the ionization dynamics of 
atomic hydrogen is mainly dictated by the Vo potential. Thus, 
in this limit, the dynamics of the system is approximately 
given by the effective Hamiltonian 

^KH = i:(y + Vo(ao,r,))+^^r^- (21) 

Note that this Hamiltonian is time-independent and accounts 
for shake-up (excitation) and shake-off (ionization) in the 
strong field limit. 

An analysis of the properties of the Vq potential term in the 
vicinity of the origin reveals that it can be neglected relative 
to the electron-electron repulsion term in the limit Oo — > °o, 
provided the two-electron wave function is localized, i.e., 
(ri) ^ Oo and (r2) <C Oq. This means that the dynamics of 
the two-electron system, in the limit of superintense fields and 
for sufficiently short pulses, ultimately reduces to that of a 
pure Coulomb explosion process effectuated by the Coulomb 
repulsion term in Eq. (l2ll . From this we conclude that the 
electron-electron interaction is in fact very important in the 
strong field limit, effectively reducing the stabilization effect. 
Returning to the laboratory frame of reference this should be 
understood in the following way: In the very strong field limit, 
the electrons effectively behave like free particles in the field, 
moving side by side with respect to the field axis. As this 
happens, the nuclear attraction may become less important 
than the mutual repulsion between the electrons, and the ion- 
ization is most likely initiated by electron-electron scattering 
events (Coulomb explosion) and not electron-nucleus colli- 
sions. This explains qualitatively why the ionization probabil- 
ities, calculated within the independent-electron model, devi- 
ate increasingly from the exact ones in the limit of stronger 
fields (c.f. Figs. [T] and |2]i. As such, the observed deviation is 
indeed a manifestation of the breakdown of the single-particle 
picture in superstrong fields. 

Notice that the effective Hamiltonian in Eq. (l2Tl i only de- 
pends indirectly on the laser frequency through the displace- 
ment ampUtude Oo, explaining why the validity of the IE 
model in Fig.|2]breaks down at approximately the same value 
of Oo independent of the laser frequency. 

Figure |3] shows the probability of excitation of helium for 
the cases considered in Fig. |2] Comparing Fig. |2]and|3]we 
observe that the decreasing ionization probability in the sta- 
bilization regime is accompanied by a corresponding increase 
in the excitation probability. Note that due to the high pho- 
ton energy, excitation is here caused by shake-up processes, 
merely demonstrating the importance of the Vq potential in the 
stabilization regime. The figure also clearly expresses the fact 
that shake-up is more important for the higher frequencies and 
that the IE model fails in describing shake-up (and shake-off) 
processes accurately, in accordance with the KH discussion 
above. 
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FIG. 3: (Color online) Excitation probabilities for the same scenario 
as in Fig. |2] The laser frequencies are o) = 10,8,6,4a.u., from top 
to bottom. The dashed lines are results from the IE model, while the 
solid lines are results from the full calculations. 



C. Analysis of angular and energy distributions 

Further insight into the strong-field behavior of helium may 
be gained by examining the energy distribution of the ejected 
electrons. In particular, imprints left by the electron interac- 
tion in the angular distribution of the outgoing electrons may 
give further clues as to its importance at the different field 
strength regimes considered here. 
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FIG. 4: (Color online) Double ionization energy distribution for an 
CO = 5 a.u. 6-cycle pulse (182 as), and for pulse amplitudes of 1 a.u. 
(upper panel), 10 a.u. (middle panel) and 20 a.u. (lower panel). 



In Fig.|4]the double ionization energy distribution is shown 
for three different pulse amplitudes, 1, 10 and, 20 a.u (from 
top to bottom). The pulse duration was fixed at 6 cycles 
(182 as), with a frequency of 5 a.u. At the lower intensity (am- 
plitude), one-photon ionization dominates as expected. As 
the intensity increases, two-photon ionization becomes promi- 
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nent, and higher-order double-electron above threshold ion- 
ization (DATI) peaks start to appear [30]. Since the one- 
photon process is highly correlated and depends critically on 
exchange of energy between the two electrons, it becomes less 
important at stronger fields, and two-photon double ionization 
takes over as the dominating channel. At the highest inten- 
sity, more structures appear in the energy spectrum, caused by 
sidebands in the pulse, and the one-photon ionization process 
has become negligible. 

The two-photon DATI component manifests itself as a 
single-peaked structure in Fig. H] in contrast to the com- 
mon double-peak structure associated with sequential ioniza- 
tion 1 30]. With the ultrashort pulse considered here, the sec- 
ond photon is absorbed before the residual ion has had time to 
relax to the ground state, but if the duration is increased to be- 
yond 20 cycles, relaxation may occur and a double-peak struc- 
ture appears (not shown here). The fact that the two peaks, 
corresponding to sequential two-photon double ionization in 
the long pulse limit, shift towards each other in the short pulse 
regime, and eventually merge into one single peak (located at 
equal energy sharing) is well known and has been studied in a 
series of papers in the weak field (perturbative) limit 
I62I1 . The results in Fig.|4]demonstrate that this feature survives 
in the superintense field regime, representing a clear departure 
from the independent-electron model Eq. (flU . 

Figure |5] shows the conditional angular distributions of 
the ejected electrons obtained at the two-photon DATI peak 
in Fig. m with equal energy sharing, and one of the elec- 
trons emitted along the polarization direction (indicated with 
an arrow in the figure). The figure clearly shows that the 
distribution has a backward-forward asymmetry even at the 
highest intensity considered, demonstrating the breakdown of 
the single-particle picture, wherein a symmetric double lobe 
(dipole) distribution would be found. The results are in ac- 
cordance with recent results obtained at weaker fields ll53l Isil 
I6OI1 . and shows that the back-to-back ejection mechanism is 
largely preserved even at very strong fields. 

From numerical studies of stabilization in atomic hydrogen, 
it is known that the stabilization phenomenon is accompanied 
by the appearance of slow electrons (T^, '63*1 . As the intensity 
is increased beyond the ionization maximum, where stabiliza- 
tion sets in, and for sufficiently short pulses, a peak structure 
near zero energy appears in the electron energy spectrum, be- 
coming increasingly dominant as the intensity becomes large. 
This may be understood from the Kramers-Henneberger anal- 
ysis above, and the importance of the Vb potential in the limit 
of short, intense pulses. In order to provide a baseline com- 
parison for the two-electron case considered here, we have 
calculated the energy distribution for ionization of He+, with 
identical pulse characteristics as those used in Fig. |4] The 
result is shown in the upper panel of Fig. |6] We note the pres- 
ence of ATI peaks, and, at the highest intensities, a slow elec- 
tron peak (SEP) near zero energy ifisl l63ll . The corresponding 



single ionization energy distribution of helium is shown in the 
lower panel of Fig. |6] and indeed, a slow electron peak is visi- 
ble. Note that the onset of slow electrons occurs at lower field 
strengths in the single ionization of helium than in He+, which 
is related to the different ionization potentials (Ip). 

Now, examining the lower panel in Fig. |4] it appears that 
slow electrons do not emerge in the double ionization process 
at this intensity. However, the ionization potential is greater 
than that for single ionization, and therefore higher intensi- 
ties are needed to reach the regime where a SEP may appear 
Since He+, with an ionization potential of Ip = 2a.u., exhibits 
an onset of slow electrons around 50 a.u., similar or possibly 
even higher field strength may be required for a SEP to appear 
in the double ionization of helium. We may, however, observe 
the onset of slow electrons by partitioning the double ioniza- 
tion energy distribution into a low- to high-energy part, and 
considering the ratio of these two, c.f. Fig. |7] In this figure, 
the ratio P{E < Ec)/P{E > Ec), Ec = 3 /2(0 - Ip, for different 
frequencies and pulse durations are shown, and in all cases, 
we observe an increase of low-energy electrons after the sta- 
bilization peak (indicated by the blue bar), however it is most 
pronounced for the shorter pulses. 



IV. CONCLUSION 

In conclusion, we have presented an in-depth analysis 
of two-electron dynamics driven by high-intensity ultrashort 
laser pulses in the xuv regime. Expanding on our earlier inves- 
tigation of correlation effects in the stabilization of helium, we 
have shown that stabilization occurs within a narrow interval 
of values of Oo, independent of frequency and pulse duration. 
This is also the point at which an independent-electron picture 
begins to break down, demonstrating the important role of the 
electron-electron interaction at high intensities. Through an 
analysis of a high-intensity limit form of the Hamiltonian, ex- 
pressed in the Kramers-Henneberger frame, this feature may 
be understood. Further indications of intense field correla- 
tion effects are found in the angular distributions, where a 
backward-forward asymmetry is found for a wide range of in- 
tensities. Finally, we have shown that slow electrons emerge 
at high intensities, as they do in one-electron systems, but at 
different intensities for single and double ionization. 
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FIG. 5: (Color online) Angular distributions for double ionization 
with equal energy sharing E\ =£2 = {2(0 — Ip)/2. The arrow in- 
dicates the fixed direction of the first electron . Solid (blue) line: 
£0 = la.u. Dashed (green) line: Eq = lOa.u. Dotted (red) line: 
£() = 20a. u. 
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FIG. 6: (Color online) Energy distributions as a function of laser field 
strength. Top panel: He+. Bottom panel: helium (single ionization). 
See text for details. 
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FIG. 7: (Color online) Ratio of slow to fast electrons for the dou- 
ble ionization process, shown for different pulse frequencies and du- 
rations, and plotted as a function of Oq. The blue bar indicates the 
region where the coiTesponding double ionization probability is max- 
imum; where stabilization sets in. 



